Constraining the helium abundance with CMB 



m 
o 
o 

-i— > 

o 
O 



(N ■ 
> ■ 

oo : 

oo . 

IT) ■ 

\o ■ 
o : 

co . 

o ■ 

Or 
i 1 

o ■ 
a: 



13 



Roberto Trotta^Q and Steen H. Hansen 2 -Q 
1 Departement de Physique Theorique, Universite de Geneve, 
24 quai Ernest Ansermet, CH-1211 Geneve 4, Switzerland 
2 Physik Institut, Winterthurerstrasse 190, 8057 Zurich, Switzerland 

We consider for the first time the ability of present-day cosmic microwave background (CMB) 
anisotropies data to determine the primordial helium mass fraction, Y p . We find that CMB data 
alone gives the confidence interval 0.160 < Y p < 0.50f (at 68% c.L). We analyse the impact on the 
baryon abundance as measured by CMB and discuss the implications for big bang nucleosynthesis. 
We identify and discuss correlations between the helium mass fraction and both the redshift of 
reionization and the spectral index. We forecast the precision of future CMB observations, and find 
that Planck alone will measure Y p with error-bars of 5%. We point out that the uncertainty in the 
determination of the helium fraction will have to be taken into account in order to correctly estimate 
the baryon density from Planck-quality CMB data. 



I. INTRODUCTION 

Our understanding of the baryon abundance has in- 
creased dramatically over the last few years. This im- 
provement comes from two independent paths, namely 
big bang nucleosynthesis (BBN) and cosmic microwave 
background radiation (CMB). Absorption features from 
high-redshift quasars allow to measure precisely the deu- 
terium abundance, D/H. Combined with BBN calcula- 
tions, this provides a reliable estimate of the baryon to 
photon ratio, rj. An independent determination of the 
baryon content of the universe from CMB anisotropies 
comes from the increasingly precise measurements of the 
acoustic peaks, which bear a characteristic signature of 
the photon-baryon fluid oscillations. The agreement be- 
tween these two completely different approaches is both 
remarkable and impressive (see details below). The time 
is therefore ripe to proceed and test the agreement be- 
tween other light elements which are also probed both 
with BBN and CMB. 

The helium abundance has been measured for many 
years from astrophysical systems. However, the error- 
bars are seemingly dominated by systematic errors which 
are hard to assess. Fortunately, the dependence of the he- 
lium mass fraction on the CMB anisotropies provides an 
independent way to measure Y p . The aim of this work 
is to present the first determination of the helium abun- 
dance from CMB alone, and to clarify the future poten- 
tial of this method. The latest CMB data are precise 
enough to allow taking this further step, and in view 
of the emerging "baryon tension" between BBN predic- 
tions from observations of different light elements |l| pos- 
sibly requires taking such a step. The advantage of using 
CMB anisotropies rather than the traditional astrophysi- 
cal measurements, is that CMB provide a clear measure- 
ment of the primordial helium fraction before it could 
be changed by any astrophysical process. On the other 
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hand the dependence of the CMB power spectrum on Y p 
is rather mild, a fact which makes it presently safe to 
fix the value of the helium mass fraction with zero un- 
ccrtainity for the purpose of CMB data analysis of other 
cosmological parameters. 

The paper is organized as follows. In section [H] we 
review the standard Big Bang Nucleosynthesis scenario. 
Section 11111 discusses the role of the helium mass frac- 
tion for cosmic microwave background anisotropies, the 
methods used and results. We discuss our forecast for 
future CMB observations in section UlI Dl and offer our 
conclusions in section ITV1 



II. BIG BANG NUCLEOSYNTHESIS 
A. The standard scenario 

The standard model of big bang nucleosynthesis 
(SBBN) has only one free parameter, namely the baryon 
to photon ratio 7710 = nfc/n 7 10 10 , which for long has been 
known to be in the range 1 — 10 |2j. Thus by observ- 
ing just one primordial light element one can predict the 
abundances of all the other light elements. 

The deuterium to hydrogen abundance, D/H, is ob- 
served by hy-a features in several quasar absorption sys- 
tems at high red-shift, D/H = 2.78±g"j£ x 10" 5 0, 
which in SBBN translates into the baryon abundance, 
rjio = 5.9±0.5. Using SBBN one thus predicts the helium 
mass fraction to be in the range 0.2470 < Y p < 0.2487. 
The dispersion in various deuterium observations is, how- 
ever, still rather large, ranging from D/H = 1.65±0.35 x 
10- 5 H to D/H = 3.98±o j? x 10" 5 0, which most prob- 
ably indicates underestimated systematic errors. 

The observed helium mass fraction comes from the 
study of extragalactic HII regions in blue compact galax- 
ies. One careful study [5| gives the value Yp = 
0.244±0.002; however, also here there is a large scat- 
ter in the various observed values, ranging from Y p = 
0.230 ± 0.003 over Y p = 0.2384 ± 0.0025 and 
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Besides the large scatter there is also the problem that 
the helium mass fraction predicted from observation of 
deuterium combined with SBBN, 0.2470 < Y p < 0.2487, 
is larger than (and seems almost in disagreement with) 
most of the observed helium abundances, which probably 
points towards underestimated systematic errors, rather 
than the need for new physics 0, [HJ • Figure is a com- 
pilation of the above measurements, and offers a direct 
comparison with the current (large) errors from CMB 
observations (presented in section III below) and with 
the potential of future CMB measurements (discussed in 
section III D). 




Astrophysical 



CMB alone 



FIG. 1: On the left (blue) we plot a few current direct astro- 
physical measurements of the helium mass fraction Y p with 
their 1 — a statistical errors, and the value inferred from deu- 
terium measurements combined with SBBN (red) (see text 
for references). On the right (green), a direct comparison 
with CMB present-day accuracy (actual data, this work; the 
errorbar extends in the range 0.16 < Yp < 0.50) and with 
its future potential (Fisher matrix forecast for Planck and a 
Cosmic Variance Limited experiment). 

The observed abundance of primordial 7 Li using the 
Spite plateau is possibly spoiled by various systematic 
effects 0, 0| . Therefore it is more appropriate to use 
the SBBN predictions together with observations to esti- 
mate the depletion factor fj — 7 L hhs /iLlppm instead of 
using 7 Li Q bs to infer the value of r\ [Tj, [TjJ H3 ■ 

The numerical predictions of standard BBN (as well as 
various non-standard scenarios) have reached a high level 
of accuracy [Hill H E HJ E3 , and the precision of 
these codes is well beyond the systematic errors discussed 
above. 



B. The role of neutrinos 

If the CMB-determined helium mass fraction turns out 
to be as high as suggested by SBBN calculations together 
with the CMB observation of VL^h 2 (as discussed above), 
this could indicate a systematic error in the present direct 
astrophysical helium observations. Alternatively, if the 
CMB could independently determine the helium value 



with sufficient precision to confirm the present helium 
observations, then this would be a smoking gun for new 
physics. In fact, one could easily imagine non-standard 
BBN scenarios which would agree with present observa- 
tions of 7710, while having a low helium mass fraction. All 
that is needed is additional non-equilibrium electron neu- 
trinos produced at the time of neutrino decoupling which 
would alter the n — p reaction. This could alter the re- 
sulting helium mass fraction while leaving the deuterium 
abundance unchanged. One such possibility would be a 
heavy sterile neutrino whose decay products include v e . 
A sterile neutrino with life-time of 1 — 5 sec and with 
decay channel v s — > v e + <f> with <p a light scalar (like a 
majoron), would leave the deuterium abundance roughly 
untouched, but can change the helium mass fraction be- 
tween AY p = -0.025 and AY p = 0.015 if the sterile neu- 
trino mass is in the range 1 — 20 MeV [2i|. A simpler 
model would be standard neutrino oscillation between a 
sterile neutrino and the electron neutrino. The lifetime 
is about 1 sec when the sterile state has mass about 10 
MeV, and the decay channel is v s — > v e + I + I (with I 
any light lepton) , and such masses and life-times are still 
unconstrained for large mixing angle |23f (related BBN 
issues are discussed in refs. |25[ l26l l27| h Such pos- 
sibilities are hard to constrain without an independent 
measurement of the helium mass fraction. 

Another much studied effect of neutrinos is the in- 
creased expansion rate of the universe if additional de- 
grees of freedom are present (for BBN), and the degener- 
acy between the total density in matter and relativistic 
particles (for CMBh This issue has recently been studied 
in detail in refs. [ToL l2*sj in view of the new WMAP re- 
sults, and we need not discuss this further here. We thus 
fix N v = 3.04 [2^. Also an electron neutrino chemical po- 
tential could potentially alter the BBN predictions |30| . 
however, with the observed neutrino oscillation param- 
eters the different neutrino chemical potentials would 
equilibrate before the onset of BBN 31], hence virtually 
excluding this possibility (see however |32j1. 



III. COSMIC MICROWAVE BACKGROUND 
A. Photon recombination and reionization 

The recent WMAP data allow one to determine with 
very high precision the epoch of photon decoupling, Zd cc , 
i.e. the epoch at which the ionized electron fraction, 
x e (z) = n e /njy, has dropped from 1 to its residual value 
of order 10~ 4 . Here n e denotes the number density of 
free electrons, while tih is the total number density of 
H atoms (both ionized and recombincd). After this mo- 
ment, photons are no longer coupled to electrons (last 
scattering), and they free stream. The redshift of de- 
coupling has been determined to be Zd oc = 108813, E3- 
which corresponds to a temperature of about 0.25 eV. 
Helium recombines earlier than hydrogen, roughly in two 
steps: around redshift z — 6000 Helll recombines to 



Hell, while Hell to Hel recombination begins around 
z < 2500 and finishes just after the start of H recom- 
bination (see e.g. (HHaEilll)- 

Denoting by nue and rif, the number densities per m 3 
of He atoms and baryons, respectively, the helium mass 
fraction is defined as Y p — Anne/nb. The baryon number 
density is related to the baryon energy density today, u>b, 
by lib — 11.3(1 + z) 3 ujb and we have njj = nb(l — Y p ). 
Usually, the ionization history is described in terms of 
x e (z) — n e /(rib(l — Y p )). However, for the purpose of 
discussing the role of Yp, it is more convenient to consider 
the quantity f e (z) = n e /rib instead, the ratio of free elec- 
trons to the total number of baryons. For brevity, we will 
call f e the free electron fraction. Once the baryon num- 
ber density has been set by fixing u>b, one can think of Y p 
as an additional parameter which controls the number of 
free electrons available in the tight coupling regime. The 
CMB power spectrum depends on the full detailed evolu- 
tion of the free electron fraction, but we can qualitatively 
describe the role of helium in four different phases of the 
ionization/recombination history (see Fig. EJ) . 

(a) Before Helll recombination all electrons are free, 
therefore f e (z > 6000) = 1 - Y p /2. 

(b) Hell progressively recombines and just before H re- 
combination begins, f e has dropped to the value 

/ e (*« 1100) = 1-Y P . 

(c) After decoupling, a residual fraction of free elec- 
trons freezes out, giving / e (30 < z < 800) = f r e es fts 
2.7 • 10- 5 y/U^/u b . 

(d) Reionization of all the H atoms gives f e (z < 20) = 
1 - Y p . 

During phase (a), the photon-baryons fluid is in the 
tight coupling regime. However the presence of ionized 
He increases diffusion damping, therefore having an im- 
pact on the damping scale in the acoustic peaks region. 
When the detailed energy levels structure of Hell is taken 
into account [37j , the transition to phase (b) is smoother 
than in the Saha equation approximation. Therefore the 
plateau with f e = 1 — Y p is not visible in Fig. [21 Be- 
fore H recombination, He atoms remain tightly coupled 
to H atoms through collisions, with the same dynam- 
ical behaviour. In particular, it is the total w& which 
determines the amount of gravitational pressure on the 
photon-baryons fluid, and which sets the acoustic peak 
enhancement /suppression. Hence we do not expect the 
value of Y p to have any influence on the boosting (sup- 
pression) of odd (even) peaks. The redshift of decoupling 
(transition between (b) and (c)) depends mildly on Y p in 
a correlated way with Ub, since the number density of 
free electrons in the tight coupling regime (just before H 
recombination) scales as n e = f e rib = n&(l — Y p ). Hence 
an increase in ujb can be compensated by allowing for a 
larger helium fraction. An analytical estimate along the 
same lines as in e.g. indicates that a 10% change in 
Yp affects z& cc by roughly 0.1%, which corresponds to 
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FIG. 2: Evolution of the number density of electrons nor- 
malized to the number density of baryons, f e — n e /nt, as a 
function of redshift for different values of the helium fraction 
Y p . The black-solid curve corresponds to the standard value 
Y p — 0.24, the red-dashed (blue-dot-dashed) to Y p — 0.36 
(Yp — 0.12). The labels (a) to (d) indicate the four different 
phases discussed in the text. 



Azdcc ~ 1. This is of the same order as the current 1-er 
errors on Zdec, obtained by fixing Y p = 0.24. 

After H recombination, the residual ionized electron 
fraction f™ s does not depend on Y p , but is inversely pro- 
portional to the total baryon density (phase (c)). As the 
CMB photons propagate, they are occasionally rescat- 
tered by the residual free electrons. The corresponding 
optical depth, t ics is given by 



tdec 



n™ B ca T dt 



(1) 



1.86-10 



-6 



(l + zf 



-dz. 



Performing the integral we can safely neglect the con- 
tribution of the cosmological constant at small redshift, 
since z^ec ^ ^A/^m- Retaining only the leading term, 
the approximated optical depth from the residual ioniza- 
tion fraction is estimated to be 



1.24-10- 6 (l + z d cc) 



3/2 



0.045, 



(2) 



independent of the cosmological parameters and of the 
helium fraction. Therefore after last scattering we do 
not expect any significant effect on CMB anisotropics 
coming from the primordial helium fraction, until the 
reionization epoch. 

Fairly little is known about the exact reionization 
mechanism and its redshift dependence (for a review see 
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e.g. |23). Observation of Gunn-Peterson troughs indicate 
that the universe was completely ionized after rcdshift 
z ~ 6, when the universe seemingly com plet ed the reion- 
ization [3^| , possibly for the second time |40j| . If tempera- 
ture information only is available, CMB anisotropics are 
sensitive only to the integrated reionized fraction, repre- 
sented by the optical depth, independent of the specific 
reionization history. However, specific signatures are im- 
printed on the E-polarization and ET-cross correlation 
power spectra by the detailed shape of the reionization 
history (for a detailed discussion, see 0, 0])- 
There are several physically motivated reionization sce- 
narios, which however cannot be clearly distinguished 
at present 0, l-Hl • In this work we use the most sim- 
ple model, the sudden reionization scenario: we assume 
that at the reionization redshift z r all the hydrogen was 
quickly reionized, thus producing a sharp rise of n e from 
its residual value to tih- More precisely, z r is the red- 
shift at which x e {z r ) = 0.5. In our treatment we neglect 
Hell reionization, for which there is evidence at a red- 
shift z ps 3 (see [42] and references therein). This effect 
is small, since the extra electron released at z « 3 would 
change the reionization optical depth by about only 1%. 
We also neglect the increase of the helium fraction due 
to non-primordial helium production, which has a negli- 
gible effect on CMB anisotropics. Those approximations 
do not affect the results at today's level of sensitivity 
of CMB data: for WMAP noise levels, even inclusion 
of the polarization spectra is not enough to distinguish 
between a sudden reionization scenario and a more com- 
plex reionization history. At the level of Planck a more 
refined model ling of the reionization mechanism will be 
necessary [44I l48j. 

In the sudden reionization scenario adopted here, the 
relation between reionization redshift and reionization 
optical depth, r r , is given by 

r r = / n e c<jTdt 

« 11.3cct t cj 6 (1 - Y p ) f ^-dz, (3) 
Jo da 

where t is physical time, 77 is conformal time and a the 
scale factor. Here again, since the number density of 
reionized electrons scales as cj&(1 — Y p ), the redshift of 
reionization is positively correlated with Y p (for fixed op- 
tical depth and baryon density). 

As a result of the physical mechanism described above, 
a 10% change in Y p has a net impact on the CMB 
power spectrum at the percent level. The impact on 
the CMB temperature and polarization power spectra 
is highlighted in Fig. |3| In the temperature panel, we 
notice that a larger helium fraction slightly suppresses 
the peaks because of diffusion damping, while it has no 
impact on large scales. Polarization is induced by the 
temperature quadrupole component at last scattering. 
When reionization occurs, there is a generation of po- 
larized power on the scale corresponding to the acoustic 
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FIG. 3: CMB temperature (top panel) and polarization (bot- 
tom panel) power spectra and percentage difference with two 
different values of the helium fraction for a standard ACDM 
model. The solid-black (dashed-blue) line corresponds to a 
10% larger (smaller) value of Y p wrt to the standard value, 
Y p — 0.24. All other parameters are fixed to the value of our 
fiducial model (Table [IJ, in particular, we have r r = 0.166. 



horizon size at the reionization redshift. This particular 
signature is called the "reionization bump" , and is clearly 
visible in the bottom panel of Fig. Oin the i ss 10 region. 
The position of the bump in multipole space scales as 
4um P oc H^. As discussed above, a change in the 
helium fraction implies a shift of the redshift of reion- 
ization for a given (fixed) optical depth, Therefore the 
value of Y p has an effect on the position of the reioniza- 
tion bump in the polarization power spectrum, but not 
on its height, which is controlled by the optical depth and 
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is proportional to t 2 . This effect is highlighted in the bot- 
tom panel: a 10% change in Y p induces roughly a 10% 
change in the position of the bump. The subsequent two 
oscillatory features for £ < 50 reflect the displacement of 
further secondary, reionization induced polarization os- 
cillations. However, since the value of polarized power is 
very low in that region, such secondary oscillations are 
very hard to detect precisely. In principle, given an accu- 
rate knowledge of the reionization history, the effect of Y p 
on the polarization bump would assist into determining 
the helium abundance. However, our ignorance of the 
reionization history prevents us from recovering useful 
information out of the measured reionization bump. The 
displacement induced by Y p is in fact degenerate with a 
partial reionization, or with other, more complex reion- 
ization mechanisms (see |3). Hence constraints on Y p 
come effectively from the damping tail in the £ > 400 
region of the temperature spectrum, which needs to be 
measured with very high accuracy. 

Other light elements like deuterium and helium-3 are 
much less abundant, and will therefore have even smaller 
effect on the CMB power spectrum, at the order of 10~ 5 . 



B. Monte Carlo analysis 

We use a modified version of the publicly available 
Markov Chain Monte Carlo package cosmomc as 
described in [EJ in order to construct Markov Chains in 
our 7 dimensional parameter space. We sample over the 
following set of cosmological parameters: the physical 
baryon and CDM densities, lu^ = Vl^ti 2 and oj c = il c h 2 , 
the cosmological constant in units of the critical density, 
f^A, the scalar spectral index and the overall normal- 
ization of the power spectrum, n s and A s (see section 
III D below for a more precise definition), the redshift 
at which the reionization fraction is a half, z r , and the 
primordial helium mass fraction, Y p . We restrict our 
analysis to flat models, therefore the Hubble parame- 
ter, h = iJo/100 km s _1 Mpc _1 , is a derived parameter, 
h = [{lj c + uJb)/(l — ^a)] 1 / 2 . We consider purely adi- 
abatic initial conditions, and we do not include gravita- 
tional waves. In the CMB analysis, we assume 3 massless 
neutrino families and no massive neutrinos. We include 
the WMAP data [H, El] (temperature and polarization) 
with the routine for computing the likelihood supplied by 
the WMAP team EH- We make use of the CBI [H and 
of the decorrelated ACBAR EH E3 band powers above 
£ = 800 to cover the small angular scale region of the 
power spectrum. 

Since Y p is a rather flat direction in parameter space 
with present-day data, we find that a much larger number 
of samples is needed in order to achieve good mixing and 
convergence of the chains in the full 7D space. We use 
M = 4 chains, each containing approximately N = 3- 10 5 
samples. The mixing diagnostic is done on the same lines 
as in [54j , by means of the Gelman and Rubin criterion 
[58| . The burn-in of the chains also takes longer than 



in the case where Y p is held fixed, and we discard 6000 
samples per chain. 

C. CMB analysis results 

Marginalizing over all other parameters, we find that 
the helium mass fraction from CMB alone is constrained 
to be Y p < 0.647 at 99% c.l. (1 tail limit), and 

0.160 < Y p < 0.501 (4) 

at 68% c.l. (2 tails). Thus, for the first time the pri- 
mordial helium mass fraction has been observed using 
the cosmic microwave background. However, present-day 
CMB data do not have sufficient resolution to discrim- 
inate between the astrophysical helium measurements, 
Y p ~ 0.244, and the deuterium guided BBN predictions, 
Y p ~ 0.248, which would require percent precision. 




FIG. 4: One-dimensional posterior likelihood distribution for 
the helium mass fraction, Y p , using CMB data only. The 
solid-black line is for all other parameters marginalized, the 
dashed-red line gives the mean likelihood. 

In Figure 0] we plot the marginalized and the mean 
likelihood of the Monte Carlo samples as a function of 
Y p . If the likelihood distribution is Gaussian, then the 
2 curves should be indistinguishable. The difference be- 
tween marginalized and mean likelihood for Y p indicates 
that the marginalized parameters are skewing the distri- 
bution, and therefore that correlations play an important 
role. Although the mean of the ID marginalized likeli- 
hood is rather high, (C(Y p )) = 0.33, the mean likelihood 
peaks in the region indicated by astrophysical measure- 
ments, Y p ~ 0.25. In view of this difference, it is im- 
portant to understand the role of correlations with other 
parameters, and we will turn to this issue now. 
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In Figure[5]we plot joint 68% and 99% confidence con- 
tours in the (u>b, Y p )-apa,ce. From the Monte Carlo sam- 
ples we obtain a small and negative correlation coeffi- 
cient between the two parameters corr(Y^, u&) = —0.14. 
Baryons and helium appear to be anticorrelated simply 
because present-day WMAP data do not map the peaks 
structure to sufficiently high I. Precise measurements 
in the small angular scale region should reveal the ex- 
pected positive correlation between the baryon and he- 
lium abundances, which is potentially important in order 
to correctly combine BBN predictions and CMB mea- 
surements of the baryon abundance. We turn to this 
question in more detail in the next section. In SBBN 
the baryon fraction and helium fraction are correlated 
along a different direction (cf. Fig. EJ. However, this 
correlation is very weak, and the SBBN relation gives 
practically a flat line. Since the two parameters are not 
independent from the CMB point of view, it is in fact not 
completely accurate to do the CMB analysis with fixed 
helium mass fraction of Y p = 0.24 to get the error-bars on 
the baryon fraction, and then re-input this baryon frac- 
tion (and error-bars) to predict the helium mass fraction 
from BBN. The most accurate procedure is to analyse 
the CMB data leaving Yp as a free parameter, thereby 
obtaining the correct (potentially larger) error-bars on 
uib upon marginalization over Y p . 



0.9- 



0.8- 
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FIG. 5: Joint 68% and 99% confidence contours in the 
(u>b, Y^)-plane from CMB data alone. The solid-blue line gives 
the SBBN prediction 0, which on this figure almost looks 
like a straight line. 

In view of the emerging baryon tension between CMB 
and BBN, it is important to check whether allowing he- 
lium as a free parameter can significantly change the 
CMB determination of the baryon density or its error. 
In order to evaluate in detail the impact of Y p on the 
error-bars for u>b, we consider the following 3 cases. 



(a) The usual case, when the helium fraction for the 
CMB analysis is assumed to be known a priori and 
is fixed to the canonical value Y p — 0.24. 

(b) A case with a weak astrophysical Gaussian prior 
on the helium fraction, which we take to be 1^ = 
0.24 ± 0.01. As discussed above, the error-bars of 
the astrophysical measurements are typically a fac- 
tor 5 tighter than this, but our prior is chosen to 
encompass the systematic spread between the dif- 
ferent observations. 

(c) The case in which we assume a uniform prior for 
Y p in the range < Y p < 1, i.e. Y p is considered as 
a totally free parameter. 

We do not find any significant change in the error- 
bars for ujb in the 3 different cases. The confidence 
intervals on u>b alone are determined to be (case (c)) 
0.0221 < uj b < 0.0245 at 68% c.l. (0.0204 < uj b < 0.0276 
at 99 % c.l.). The standard deviation of LOb as esti- 
mated from the Monte Carlo samples is found to be 
(7/, = 1.3 • 10~ 3 . This is in complete agreement with the 
error-bars on u>b obtained by the WMAP team for the 
standard ACDM case |33| . We conclude that at the level 
of precision of present-day CMB data, it is still safe to 
treat the baryon abundance and the helium mass fraction 
as independent parameters. This result is non-trivial, 
since the fact that the damping tail is not yet precisely 
measured above the second peak would a priori suggest 
that degeneracies between Y p , uj b and n s could potentially 
play a role once the assumption of zero uncertainity on 
Y p is relaxed. The impact of Y p is small enough, and the 
error-bars on u>b large enough that a uniform prior on Y p 
can still be accomodated within the uncertainity in the 
baryon abundance obtained for case (a). However, the 
Y p — LUb correlation will have to be taken into account to 
correctly analyse future CMB data, with a quality such 
as Planck. We discuss this potential in the next section. 

We observe the expected correlation between the red- 
shift of reionization and the helium fraction (Fig. EJ, 
which is discussed above. The correlation coefficient be- 
tween the two parameters is found to be rather large and 
positive, corr(Yp,z r ) = 0.40. This correlation produces 
a noticeable change in the marginalized ID-likelihood 
distribution for z r as we go from case (a) to case (c). 
Marginalization over the additional degree of freedom 
given by Y p broadens considerably the error-bars on z r . 
In fact, the 68% confidence interval for z r increases by 
roughly 20% (and shifts to somewhat higher values), from 
10.2-20.9 (case (a)) to 10.6-23.3 (case (c)). Case (b) ex- 
hibits similar error-bars as case (a). On the other hand, 
the determination of the reionization optical depth is not 
affected by the inclusion of helium as a free parameter, 
giving in all cases 0.08 < r r < 0.23. Correspondingly, the 
correlation is less significant, corr(Yp,T r ) = —0.11. We 
therefore conclude that the differences in the determina- 
tion of z r are due only to the variation of the amount of 
electrons available for reionization as Y p is changed. 
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FIG. 6: Joint 68% and 99% confidence contours in the 
(Y p , z r )-plane (upper panel) and in the corresponding (Y p ,t t )- 
plane (bottom panel) from CMB data alone. In the upper 
panel, the solid-red line is the relation z r (Y p ) from eq. l|3J, 
obtained by fixing the reionization optical depth to the value 
r r = 0.166, while the other parameters are the ones of our 
fiducial ACDM model. Although clearly the exact shape of 
z r (Y p ) depends on the particular choice of cosmology, it is ap- 
parent that the Y p — z r degeneracy is along this direction. The 
correlation between Y p — r r is almost negligible with present- 
day data (bottom panel). 



Leaving Y p as a free parameter also has an impact on 
the relation between u>b and the scalar spectral index, 
n s . The extra power suppression on small scales which 
is produced by a larger Y p can be compensated by a blue 
spectral index (see Fig.0). 
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FIG. 7: Scatter plot in the u>b — n s plane, with the value 
of Y p rendered following the colour scale. Green corresponds 
roughly to the SBBN preferred value. 



described in terms of a few experimental parameters. 
The covariance matrix C is then given by the inverse 
of the Fisher matrix, C = F -1 . It is then straightfor- 
ward to evaluate the expected 1-er error on parameter 
i, which is given by yjcii (all other marginalized). The 
main advantage of the FMA is that it gives reliable and 
accurate predictions (including information on the ex- 
pected degeneracies) with minimal computational effort. 
For further details on the Fisher Matrix formalism, see 

e.g. m m mm him m m. 



1. Parameters set 



D. Potential of future CMB observations 

In order to estimate the precision with which future 
satellite CMB measurements will be able to constrain the 
helium mass fraction we perform a Fisher matrix anal- 
ysis (FMA). This technique approximates the likelihood 
function with a Gaussian distribution around a fiducial 
model, which is assumed to be the best fit model. The 
Fisher information matrix F gives the second order ex- 
pansion of the likelihood around its peak, and it is com- 
puted from the derivatives of the power spectrum with 
respect to the cosmological parameters. The expected 
performance of the experiment can be modelled with a 
noise contribution to the likelihood function, which is 



In order to obtain a reliable prediction, it is extremely 
important to choose a parameter set wrt which the de- 
pendence of the CMB power spectrum is as linear and 
uncorrelated as possible. This issue has been discussed 
exhaustively in Ref. |67l |. where the authors introduce a 
set of "physical parameters" which satisfies the above re- 
quirements. In the present work we retain most of the 
physical parameters defined in ref. |67| : the ratio between 
the sound horizon at decoupling and the angular diam- 
eter distance A, the baryon density B = ri&/i 2 , the en- 
ergy density in the cosmological constant V = fl\h 2 , the 
matter-radiation density at decoupling 1Z and M. , which 
is mainly a function of the matter and radiation con- 
tent. We adopt a slightly different choice for the physical 
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parameter describing reionization. For adiabatic pertur- 
bations, the initial power spectrum of the gauge invariant 
curvature perturbation ( is written as 

s) 

(and we do not take running into account). The quan- 
tity £ corresponds to the intrinsic curvature perturba- 
tion on comoving hypersurfaces, and at the end of in- 
flation is related to the gravitational potential perturba- 
tion, ^, by ( — |\I> (see e.g. [68| for more details). We 
take the pivot-scale fco to be ko — 0.05 Mpc . If r r 
denotes the optical depth to reionization, then defining 
T = A s exp(— 2r) is a good way to take into account 
the degeneracy between the optical depth and normal- 
ization. Our parameter set contains then the six above 
physical parameters (A,B,V,TZ,M,T), the power spec- 
trum normalization A s , the scalar spectral index n s and 
the helium fraction Y p . 

The choice of the physical parameter set makes it easy 
to implement in the FMA interesting theoretical priors. 
For instance, we are interested in imposing flatness in our 
forecast, in order to be able to directly compare present- 
day accuracy on Y p with the potential of Planck and 
CVL. The prior on the curvature of the universe is im- 
posed in the FMA by fixing the value of the parameter A 
to the one of the fiducial model. In fact, the parameter 
A is a generalization of the shift parameter, which de- 
scribes the sideway shift of the acoustic peak structure of 
the CMB power spectrum as a function of the geometry 
of the universe and its content in matter, radiation and 
cosmological constant. Although imposing A = const is 
not the same as having curvature=constant over the full 
range of cosmological parameters, for the purpose of eval- 
uating derivatives the two conditions reduce to the same. 
The fact that our fiducial model is actually slightly open 
(see below), does not make any substantial difference in 
the results, apart from reducing the numerical inaccura- 
cies which would arise had we computed the derivatives 
around an exactly flat model. We can also easily impose 
a prior knowledge of the helium fraction, by fixing the 
value of Y p , as it is usually done in present CMB analy- 
sis, and investigate how this modifies the expected error 
on the the baryon density. 

2. Accuracy issues 

We numerically compute double sided derivative of the 
power spectrum around the fiducial model with cosmo- 
logical parameters given in Table [Q We find it neces- 
sary to increase the accuracy of CAMB by a factor of 
3 in each of the "accuracy boost" values. As a fiducial 
model, we use the best fit model to the WMAP data 
for the standard ACDM scenario, as given in Table 1 
of ref. 33] . However, in order to avoid numerical inac- 
curacies which arise when differentiating around a flat 



TABLE I: Cosmological parameters for the fiducial ACDM 
model around which the FMA is performed. We choose a 
slightly open model to avoid numerical inaccuracies in the 
derivatives. 



Baryons 


n b 


0.046 


Matter 




0.270 


Dark Energy 


Ha 


0.720 


Radiation 


^rad 


7.95 • 10" 3 


Massless v families 


N v 


3.04 


Total density 


"tot 


0.990 


Hubble constant 


h 


0.72 


Optical depth 


Tr 


0.166 


Spectral index 


n 3 


0.99 


Normalization 


A s 


2 ■ 10 -9 



model, we reduce slightly the value of by imposing 
an open universe, f2 tot = 0.99. We perform the FMA 
for the expected capabilities of Planck's High Frequency 
Instrument (HFI) and for an ideal CMB measurement 
which would be cosmic variance limited (CVL) both in 
temperature and in E-polarization (and we do not con- 
sider the B-polarization spectrum), and therefore rep- 
resents the best possible parameter measurement from 
CMB anisotropies alone. The complicated issues coming 
from foreground removals, point source subtractions, etc. 
are assumed to be already (roughly) taken into account in 
the experimental parameters for the experiment. Those 
are the effective percentual sky coverage / 8 k y , the num- 
ber of channels, the sensitivity of each channel aJ ,E for 
temperature (T) and E-polarization (E) in [iK and the 
angular resolution Q^ ,E (in arcmin). For Planck HFI, 
we take the 3 channels with frequencies 100, 143 and 
217 GHz, with respectively cr^ =1 2 3 = 5.4,6.0,13.1 and 
a c=2 3 = 11-4, 26.7 and we have / s k y = 0.85 [6!j Since the 
CVL is an ideal experiment, we put its noise to zero and 
assume perfect foregrounds removal, so that / s k y = 1. In 
order to test the accuracy of our predictions and compare 
present-day results with the forecasts, we also perform 
an FMA with WMAP first year parameters, obtaining 
excellent agreement between the FMA results and the 
error-bars from actual data. For the purpose of compari- 
son, we include forecasts for the full WMAP 4 years mis- 
sion, which will also measure E-polarization and reduce 
present-day errors on the temperature spectrum by a fac- 
tor of 2. We limit the range of multipoles to £ < 2000, be- 
cause at smaller angular scales non-primary anisotropies 
begin to dominate (Sunyaev-Zeldovich effect). The au- 
thors of ref. discuss the issue of numerical precision 
of 3 different CMB codes and conclude that they are ac- 
curate to within 0.1%. While this is encouraging, it is 
not of direct relevance to this work, since what matters 
in the computation of derivatives is not much the ab- 
solute precision of the spectra, but rather their relative 
accuracy. 
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TABLE II: Fisher Matrix forecasts and comparison with 
present-day results, for different priors and using different 
combinations of temperature and polarization CMB spectra. 
Errors are in percent wrt the values of the fiducial model, 
Y p — 0.24 and uib = 0.0238 ( 1-<t c.l. all other marginalized). 



Temperature, TE-cross, E-polarization 




No priors 


Flatness 


Flatness and 




Y p w b 


AYp Au b 


Y p = 0.24 


WMAP 4yrs a 


~ 50 2.92 


~ 40 2.86 


2.86 


Planck 


7.60 1.31 


4.96 1.26 


0.70 


CVL 


2.59 0.34 


1.52 0.32 


0.13 


Temperature + TE-cross 


WMAP 1st yr b 


N/A N/A 


71.25 5.04 


5.04 


WMAP 4yrs a 


~ 75 4.10 


~ 60 3.94 


3.94 


Planck 


8.91 1.74 


6.60 1.63 


0.74 


CVL 


5.18 0.55 


2.84 0.55 


0.19 



a FMA forecast, 4 years mission including E-polarization. 

b Actual WMAP data and other CMB experiments, this work. 



3. Forecasts and discussion 

Table HTl summarizes our forecasts for the future mea- 
surements and compares them with the results obtained 
from WMAP actual data. 

We notice that when the WMAP full 4 years data will 
be available (including E-polarization), the error on the 
baryon density is expected to decrease by a factor of 
2 to 2.86%, compared to today's 5.04% (assuming flat- 
ness). Nevertheless, inclusion of Y p as a free parameter 
will still have no effect on the determination of u>f, for 
WMAP, i.e. Y p will remain an essentially flat direction 
when marginalized over. While the determination of the 
helium fraction will improve, the FMA cannot reliably 
assess quantitatively how much, since for such large er- 
rors the likelihood distribution is not Gaussian and the 
quadratic approximation breaks down. In the table we 
therefore give the FMA estimation as an indication, with 
the caveat that the Fisher approximation is likely to be 
inaccurate for the real errors on Y p from WMAP's 4 years 
data. 

It is interesting that for Planck the effect of the he- 
lium fraction can no longer be neglected. Inclusion of 
the helium fraction increases the error on tub by roughly 
80%, from 0.70% to 1.26%. The correlation between the 
two parameters will have to be taken into account, as is 
evident from Fig. EI The expected correlation coefficient 
is corr(y p ,cJ 6 ) = 0.84 (0.91) for Planck (for CVL). The 
expected 1 — a error on Y p is about 5% for Planck, or 
AY P ~ 0.01. This is of the same order as the spread in 
current astrophysical measurements. We conclude that 
in Planck-accuracy data analysis it will be necessary to 
include the uncertainty in the determination of the he- 
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FIG. 8: FMA forecast for the expected errors from WMAP 
4 years mission (dotted-black) , Planck (dashed-red) and a 
CVL experiment (solid-green). The ellipses encompass 1-cr 
and 3-cr joint confidence regions for u>b — Y p (all other param- 
eters marginalized). The axis values give the error in wrt the 
fiducial model values. This forecast is for the full CMB infor- 
mation (Temperature, TE-cross, E-polarization) and assumes 
flatness. 



hum mass fraction, at least in the form of a Gaussian 
prior over Y p of the type we used in the CMB data anal- 
ysis presented above. 

Finally, measuring CMB temperature and polarization 
with cosmic variance accuracy would allow to constrain 
Y p to within 1.5%, or AY p ~ 0.0036 (assuming flatness). 
Such an ideal measurement would be able to discrimi- 
nate between the BBN-guided, deuterium based helium 
value and the current lowest direct helium observations 
(cf. Fig. HI). 

Our forecasts for the uncertainty in the Helium mass 
fraction from future observations are in excellent agree- 
ment with the findings of Ref. There, the standard 
deviation on Y p for Planck is estimated to be AY p = 
0.012. The authors of Ref. [fj also consider an experi- 
ment (CMBPol) with characteristics similar to our CVL, 
for which they forecast AY p = 0.0039, again in close 
agreement with our result. An earlier work |72| found 
for Planck (temperature and polarization) AY p = 0.013, 
in satisfactory concordance with our result. It should be 
noticed that the forecast reported for MAP in Table 2 of 
Ref. |x3 1 namely AY p = 0.02, is nothing but the Gaus- 
sian prior Y p = 0.24 ± 0.02 which was assumed in their 
analysis. 

The main source of improvement for the determina- 
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tion of Y p will be the better sampling of the temperature 
damping tail provided by Planck and the CVL. Polar- 
ization measurements have mainly the effect of reducing 
the errors on other parameters. In fact, we have checked 
that excluding from our FMA the 2 < i < 50 region of 
the E-polarization and ET-correlation spectra changes 
the forecast precision on Y p less than about 10-15% for 
Planck and less than a few percent for CVL. This sup- 
ports the conclusion that the low-£ reionization bump is 
not very useful in measuring the helium abundance, be- 
cause of the degeneracy with z r . 



5% (or AY p ~ 0.01), which however will only allow a 
marginal discrimination between different astrophysical 
measurements. Nevertheless, we point out that the un- 
certainty of the helium fraction will have to be taken 
into account in order to correctly estimate the errors on 
the baryon density from Planck. To determine if the 
emerging baryon tension (from BBN) is related to un- 
derestimated systematic error-bars or whether it is an 
indication of new physics, CMB observation will have to 
be pushed very close to the cosmic variance limit in both 
temperature and polarization. 



IV. CONCLUSIONS 

We have analysed the ability of CMB observations to 
determine the helium mass fraction, Y p . We find that 
present data only allow a marginal detection, 0.160 < 
Y p < 0.501 at 68% c.L. This determination is com- 
pletely independent from the usual astrophysical obser- 
vations and uses CMB data only. We discuss degenera- 
cies between Y p and other cosmological parameters, most 
notably the baryon abundance, the redshift and optical 
depth of reionization and the spectral index. We con- 
clude that present-day CMB data accuracy does not re- 
quire the inclusion of Y p as a free parameter. We find that 
Planck will determine the helium mass fraction within 
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